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ABSTRACT 

We analyze the latent heat and surface tension at the SU(3) deconfinement phase 
transition with tree level and tadpole improved Symanzik actions on lattices with 
temporal extent Nr = 3 and 4 and spatial extent N^r/Nr = 4, 6 and 8. In comparison 
to the standard Wilson action we do find a drastic reduction of cut-off effects already 
with tree level improved actions. On lattices with temporal extent = 4 results for 
the surface tension and latent heat obtained with a tree level improved action agree 
well with those obtained with a tadpole improved action. A comparison with A^^ = 
3 calculations, however, shows that results obtained with tadpole action remain 
unaffected by cut-off effects even on this coarse lattice, while the tree level action 
becomes sensitive to the cut-off. 



For the surface tension and latent heat we find (Ji/Tj^ = 0.0155 (16) and Ae/T^ = 
1.40 (9), respectively. 



1 Introduction 



Tree level improved Symanzik actions [|I| have been shown to yield a large reduction 
of finite cut-off effects on bulk thermodynamic observables even at temperatures 
close to Tc ||2|, ||. It has been found that the large cut-off dependence observed in 
calculations with the standard Wilson action 0] to a large extent does result from 
cut-off dependent corrections to the high temperature ideal gas limit. They, there- 
fore, can be drastically reduced already with tree level improved Symanzik actions 
0, ^ or fixed point actions 0, which are constructed in order to reduce or elimi- 
nate the cut-off dependence in the weak coupling / high temperature limit. Further 
cut-off dependent corrections of 0{a'^g'^) do seem to be small for thermodynamic 
observables in the high temperature plasma phase. This can be seen, for instance, 
from a calculation of the difference of energy density and three times the pressure. 
The leading ideal gas contributions cancel in this quantity. It thus is sensitive only 
to 0{g'^) and higher order corrections. Calculations of this quantity performed with 
the tree level improved Symanzik action on lattices with small temporal extent do 
agree with the continuum extrapolation obtained from simulations with the Wilson 
action . Bulk thermodynamics of SU (N) gauge theories can thus be studied with 
tree level improved actions showing little cut-off distortion already on lattices with 
temporal extent = 4. 

The success of improved actions for the calculation of bulk thermodynamics even 
at temperatures close to Tc naturally leads to the question whether these actions 
also do lead to an improvement at Tc. This is, of course, a highly non-perturbative 
regime. However, observables like the latent heat and the surface tension, which 
characterize the discontinuities at the first order deconfinement phase transition 
in a SU{3) gauge theory, do depend on properties of the low as well as the high 
temperature phase. As the latter is largely controlled by high momentum modes it 
may be expected that some improvement does result already from a perturbatively 
improved description of the high momentum modes in the deconfined phase. We 
will show here that this is indeed the case. In addition we will address the question 
to what extent a tadpole improved Symanzik action ||^, |^ does lead to a further 
reduction of the cut-off dependence and would allow to perform calculations on 
even coarser lattices. 

The latent heat and the surface tension at the SU(3) deconfinement transition 
have been studied with the standard Wilson action on lattices up to temporal ex- 
tent A^T- = 6 0. A strong cut-off dependence has been found when comparing 
calculations for = 4 and 6 which is compatible in magnitude with the cut-off 
dependence visible in calculations of the energy density or pressure in the plasma 
phase. On the basis of these results an extrapolation to the continuum limit for 
these observables is not yet possible and calculations on larger temporal lattices do 
not seem to be feasible with the standard Wilson action given the present computa- 
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tional possibilities. One thus has to aim at a reduction of the cut-off dependence on 
lattices with small temporal extent in order to get reliable results for the continuum 
values of the latent heat and surface tension at Tc. We will present here results 
from calculations with a tree level as well as a tadpole improved Symanzik action 
on lattices of temporal extent Nr = 3 and 4 and spatial extent Nfj/Nr ranging from 
4 to 8. 



2 Critical Couplings 

We use the (1,2) Symanzik action which is defined by adding a planar (1,2) Wil- 
son loop to the standard Wilson action defined in terms of a (1,1) Wilson loop 
(plaquette), 

s(^) - Y.^(^Kii:'=)-^Ki(^)) . (2.1) 

Here W^'l denotes a symmetrized combination of A; x / Wilson loops in the (/x, v)- 
plane of the lattice, 

w'.fA^) = 1 - ^ (Re Trv;(5v;«,,,,v;(5t^\/if + {k^ /)) , (2.2) 

with V,^^^ = Yi^Zo Ux+j[i,n and x = (ni, n2, n^, n^j denoting the sites on an asymmet- 
ric lattice of size N^N^. We consider the tree level improved action, corresponding 
to Uo = 1, as well as a tadpole improved action |^ where 

-0 = ^(E(l-O-))) (2.3) 

is the self-consistently determined plaquette expectation value. We have deter- 
mined this tadpole improvement factor at several values of the gauge coupling (3 
and then used a spline interpolation in order to define Uq at intermediate values of 
j3. This allows to give an unambiguous definition of derivatives of the tadpole term 
with respect to the gauge coupling which is needed in thermodynamic observables. 
Moreover, it allows us to use the Ferrenberg-Swendsen reweighting technique in the 
vicinity of Tc. We note that the tadpole improved action does depend on the gauge 
coupling [3 through uq. This does lead to some modifications of the relations usually 
explored for thermodynamic observables, which are defined through derivatives of 
the partition function, Z = J VU exp{—pS{P)), with respect to /5. For the purpose 
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of the calculations presented here this is of relevance for the analysis of the latent 
heat. 

Most of our calculations have been performed close to the critical coupling on lat- 
tices with temporal extent = 3 and 4. On the smaller lattices we have performed 
simulations at several nearby /3-values. Typically 20.000-40.000 configurations have 
been generated at each /5- value using an overrelaxed heat bath algorithm (1 itera- 
tion = 4 overrelaxation + 1 heat bath step). The critical couplings have then been 
determined from the location of peaks in the Polyakov-loop susceptibilities, 

X=W)-{\L\r , (2.4) 



where L denotes the Polyakov-loop 

^ = 4ETr n %,^4),4 . (2.5) 

cr X X4 = l 

We have used an interpolation based on the Ferrenberg-Swendsen reweighting method 
to determine the location of the maximum in x- In Table |l| we summarize our results 
for the critical couplings. In the case of the Symanzik improved tree level action our 
results do agree well with those of earlier calculations For A^^ = 4 we have used 
results from different spatial lattice sizes to perform an infinite volume extrapolation 
of the critical couplings using the ansatz 

Pc{Nr) = Pcioo) - h{K/N„f (2.6) 



which is valid for first order phase transitions. The coefficient h characterizing the 
infrared sensitivity of the thermodynamics at has been found to be similar for 
both improved actions. We find h = 0.101 (34) and 0.068 (45) for the tree level 
and tadpole improved actions, respectively. This also agrees well with the volume 



dependence observed in the case of the standard Wilson action P, 11 



3 Surface Tension 



The SU{3) deconfinement transition is known to be of first order. At the critical 
temperature the deconfined and confined phases can coexist in a mixed phase. Both 
coexisting phases are separated by an interface. The additional energy needed to 
build up such interfaces does, however, make these configurations more unlikely 
than pure states of one or the other phase. This is reflected in the double-peaked 
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tree level 


tadpole 








# iterations 






^ iterations 


12^ X 3 


3.9079 (6) 


5 


80000 


4.1868 (4) 


6 


85500 


16^ X 4 


4.0715 (4) 


8 


157500 


4.3512 (4) 


7 


133400 


24^ X 4 


4.0722 (4) 


3 


75200 


4.3519 (5) 


1 


38500 


32^ X 4 


4.0729 (3) 


1 


46800 


4.3522 (4) 


1 


59000 


(oo)^ X 4 


4.0730 (3) 






4.3523 (4) 







Table 1: Critical couplings for the tree level and tadpole improved actions. Errors 
are obtained from a jack-knife analysis. Columns four and seven give the combined 
statistics from different /3-values entering the Ferrenberg-Swendsen analysis of the 
susceptibilities, surface tension and latent heat. The number of /5- values entering 
this analysis are given in columns three and six. 

probability distribution functions of observables at the critical point. While the two 
maxima are proportional to the probability of finding the system in one of the two 
pure states the minimum inbetween is related to the probability of a state containing 
two phases separated by an interface. 

To be specific we consider the probability distribution of the absolute value of 
the Polyakov loop, \L\, and follow the analysis presented in Ref. |^ to extract the 
surface tension. The probability distribution at the minimum is proportional to 

P{\L\) ~ exp{-[f,V, + hV2 + 2aiA]/T) (3.1) 

where denotes the free energy in the phase and Vi is the volume occupied by 
that phase. The last term denotes the extra free energy needed to create a surface 
to separate both phases. Here aj denotes the surface free energy and A is the area 
of the interface. An additional factor of two does appear here because we consider 
finite volumes with periodic boundary conditions and one always has to create two 
interfaces in this case. At the free energy density in both phases is identical, 
/i = /2- In our simulations we have observed frequent flips between the confined 
and deconfined phase on the 12^, 16^ and 24^ lattices which suggests that we have 
properly sampled both phases. On the 32^ lattice the data samples include, however, 
only 4-5 flips at each /5- value. 

In order to determine the maxima and the minimum of P{\L\) the region in 
the vicinity of the extrema has been fitted with third order polynomials A 
Ferrenberg-Swendsen reweighting has been used to shift the /3- value to a point where 
the two maxima, Pmax,i and Pmax,2-, in the distribution function of the Polyakov-loop 
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have equal height. A simple estimate of the surface tension is then given by 




0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 




0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 

Figure 1: Probabihty distributions of \L\ at the critical couplings. 




where Pmax = Pmax,i = Pmax,2- In Order to obtain a reliable error estimate we have 
performed a jack-knife analysis at this /9-value. As the relative maxima are then no 
longer of equal height in the course of the analysis, we define the surface tension 
through the jack-knife average of 
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where 71 and 72 characterize the relative weight of the confined and deconfined 
phases at jSc and are fixed by demanding that the Polyakov loop expectation value is 
given by the weighted sum of the values Li and L2 at Pmax,i and Pmax,2, respectively, 
i.e. 

= 7iLi + 72L2 with 7i + 72 = 1 • (3.4) 

Of course, we do find that the surface tension defined this way agrees within errors 
with the global value extracted from Eq. |3.2| . For the relative weight 71 we find in 
all cases values close to 0.4 which also is similar to the case of the standard Wilson 
action^ 0. 

The results for (Ti/T^ are given in Table |^ for the improved actions analyzed by 
us on lattices with temporal extent = 3 and 4. Also given there are the results 
from P] obtained with the standard Wilson action on iV,- = 4 and Nr = 6 lattices. 











volume 


tree level 


tadpole 


Wilson 


3 


12^ 


0.0234 (24) 


0.0158 (11) 




4 


122 X 24 






0.0241 (27) 


4 


242 X 36 






0.0300 (16) 


4 


16^ 


0.0148 (16) 


0.0147 (14) 




4 


243 


0.0136 (25) 


0.0119 (21) 




4 


32^ 


0.0116 (23) 


0.0125 (17) 




4 


00 


0.0152 (26) 


0.0152 (20) 


0.0295 (21) 


6 


20^ 






0.0123 (28) 


6 


243 






0.0143 (22) 


6 


36^ X 48 






0.0164 (26) 


6 


00 






0.0218 (33) 



Table 2: Surface tension for the improved actions on several finite spatial lattices 
with temporal extent A^^ = 4 and the Wilson action for = 4 and 6. The latter 
are taken from Ref. [9]. The second column gives the spatial lattice size. 

In order to extract the infinite volume result for the surface tension one has to 
take into account finite volume corrections, which result from the finite width of 

^ There is a misprint in the labeHng of the last column of Table I in Ref. . This column does 
give 72 rather than 71 . 
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the (gaussian) peaks in the pure phases, the contribution of zero modes resulting 
from the translational invariance as well as from contributions of fluctuations of the 
interface. Taking these into account, the interface tension should be extrapolated 
to infinite volume using the ansatz given in Ref. 0, 



T3 



c+-lnN^ 



(3.5) 



The result of such an extrapolation is also given in Table ^ for the A^^ = 4 calcula- 
tions. Even without performing the infinite volume extrapolations, it is, however, 
obvious already from Table ^ that the surface tension is reduced in the improved 
action calculations relative to that of the Wilson action with the same temporal 
extent. 

We note that the tree level improved action still shows some cut-off dependence 
when comparing the A^^ = 3 and 4 results while the results for the tadpole improved 
action coincide within errors. Combining both results for the tadpole improved 
action we find as an estimate for the surface tension 



^ = 0.0155 ±0.0016 



(3.6) 



It is interesting to note that this agrees well with an extrapolation in 
of the results obtained with the Wilson action for = 4 and 6. 



(aTc 



4 The Latent Heat 



The latent heat is calculated from the discontinuity in the energy density, e, or more 
conveniently directly from the discontinuity in (e — 3P). The latter is obtained from 
the discontinuity in the various Wilson loops entering the definition of the improved 
actions. 



with S = S ~ dS/dp and adp/da denoting the SU(3) /^-function. 

The difference of the expectation values at Pc is obtained by calculating expec- 
tation values separately in the two coexisting phases at the critical coupling Pc- 
This does require large spatial lattices in order to separate well the two coexisting 
phases. We therefore have performed the analysis of the latent heat only for the 
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32'^ X 4 lattices. Following the approach used in Ref. [Q we have analyzed the time 
histories of the Polyakov loop values and introduced a cut to classify configurations 
belonging to either of the two phases. We then have performed averages in the two 
phases separately and left out configurations belonging to the transition region. We 
have checked that the expectation values calculated in both phases are within errors 
insensitive to the precise choice of the cut and the width of transition regions. In 

0.09 
0.08 
0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 


-0.004 -0.002 0.000 0.002 0.004 

Figure 2: Probability distributions of S at the critical couplings. 

Figure ^ we show the distribution of S for the improved actions as well as the Wil- 
son action on lattices with temporal extent = 4. It is obvious from this Figure 
that the discontinuity is drastically reduced for the improved actions relative to the 
Wilson action. 

In order to extract the latent heat we still need the /3-function entering the 
definition of Ae/T^ in Eq. [4.1| . We have determined this through the calculation of 
the string tension in units of the cut-off following the approach discussed in Ref. . 
These calculations have been performed on lattices of size 16^ |^. We also have 
reanalyzed the results for the latent heat presented in 0] using the newly determined 
/9-function of Further details of our determination of the temperature scale will 
be presented in Ref. 0. 

The resulting values for the latent heat are summarized in Table ^. There we 
also give results obtained with the perturbative /^-function, adp/da = 33/47r^. The 
ratio of the fourth and fifth column in Table ^ thus reflects the deviations of the non- 
perturbative determination of the /3-functions from the perturbative result. We note 
that these ratios are, in fact, quite similar, which shows that the improved actions do 
not lead to an improved asymptotic scaling behaviour. The results summarized in 
Table |^ do, however, show that the improved actions do reduce drastically the cut- 
off dependence of the latent heat. Results obtained with the tree level and tadpole 
improved actions for A',- = 4 agree with each other. They are about 30% smaller 
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action 


K 


Nr 




\ / pert 


standard Wilson 


242 X 36 


4 


2.27 (5) 


4.06 (8) 




36^ X 48 


6 


1.53 (4) 


2.39 (6) 


(1,2) Symanzik (tree level) 


32^ 


4 


1.57 (12) 


2.28 (8) 


(1,2) Symanzik (tadpole) 


32^ 


4 


1.40 (9) 


1.94 (8) 



Table 3: Latent heat for the three improved actions and the Wilson action. Re- 
sults for the Wilson action are taken from [8] using the non-perturbative /3-function 
calculated in Ref. [4]. 

than in the case of the Wilson action and are compatible with the A^,- = 6 results 
for this action. As a best estimate for the latent heat we obtain 

1^ = 1.40 ±0.09 . (4.2) 



5 Conclusion 



We have analyzed the surface tension and the latent heat at the first order decon- 
finement phase transition of the SU(3) gauge theory using tree level and tadpole 
improved actions. We do find that these actions drastically reduce the cut-off de- 
pendence observed previously in calculations of these quantities with the standard 
Wilson action. Results obtained with a Symanzik improved tree level action agree 
within errors with those obtained with a tadpole improved action for A^^ = 4 lat- 
tices. The results obtained with the tree level action do, however, show a clear 
cut-off dependence when comparing calculations on the A^,- = 3 and 4 lattices while 
results obtained with the tadpole improved action remain unaffected by the change 
in lattice cut-off. We thus expect that our results obtained on the A^^ = 4 lattice 
with the tadpole improved action do provide a good estimate for the continuum 
limit result of discontinuities in thermodynamic observables at Tc. 

Over the small range of couplings we have explored in this study the tadpole 
improvement factor varies only slightly, 0.86 < uq{i3) < 0.88. Although we do 
take into account this variation in our simulations we note that we do effectively 
work with an over- improved Symanzik action where the weight of the contribution 
of (1 X 2)-loops relative to the (1 x l)-plaquette term is increased from 1/10 by 
about 30%. Testing the role of the /3-dependence of the tadpole improvement factor 
does require an analysis of the thermodynamics in a larger temperature interval. 
This will be presented elsewhere [§]. It also would be interesting to investigate the 
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sensitivity of the analysis presented here to even larger changes in the weight of the 
(1 X 2)-loop contribution. In the case of the renormalization group improved action 
p^ , [13|, for instance, the relative weight is as large as 0.1815. 
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